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Abstract 

In present work, we theoretically study the electron wave's focusing phenomenon in a single 
layered graphene pn junction(PNJ) and obtain the electric current density distribution of graphene 
PNJ, which is in good agreement with the qualitative result in previous numerical calculations 
[Science, 315, 1252 (2007)]. In addition, we find that for symmetric PNJ, 1/4 of total electric 
current radiated from source electrode can be collected by drain electrode. Furthermore, this ratio 
reduces to 3/16 in a symmetric graphene npn junction. Our results obtained by present analytical 
method provide a general design rule for electric lens based on negative refractory index systems. 

PACS numbers: 81.05.Uw, 42.25.Fx 
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The theory of negative refraction has been firstly analyzed by Veselago in 1968^. It was 
suggested that "perfect" optical lens can be made based on negative refraction, which can 
focus light into a fine point^. Inspired by the analogy between ballistic electron transport 
in graphene and light rays in dielectric medium, Cheianov et al. proposed that an electron 
equivalent of negative refractory can be realized in a mono-layered graphene^. By fine- 
tuning the densities of charge carrier on both sides of graphene pn junction (PNJ) to be 
equal values, the electron flow radiated from a point-like electric current source can focus 
exactly on the point in the other side of PNJ due to the negative refractory of electron wave 
crossing from conduction band to valence band. The long electron mean free path, ballistic 
electronic transport, and high current density of graphene make graphene a good candidate 
for new devices based on electric lens effect.-"— 

In addition to PNJ, the electric lens of graphene npn junction (NPNJ) is quite attractive 
for its potential applications in nanoelectronics^. Electron Veselago lens of NPNJ can be 
made by a graphene with two gates, one bottom gate with positive voltage is used to provide 
N region in a mo no- layer graphene sheet, while one top gate with negative voltage ensures 
the center of the graphene sheet is P region. By carefully controlling the top gate voltage, 
the charge density of P region which is between two identical N regions can be adjusted to 
be ph = Pe or pe 7^ ph, therefore affect focusing. A graphene transistor can be achieved in 
NPNJ with a rectangle-shaped top gate electrode, and a beam splitter can be made using 
prism-shaped top gate^, which can be used to fabricate logical gates and interconnects for 
a next generation electronics. The principle of graphene NPNJ for point-like electric source 
is similar to that of PNJ, but has not been studied quantitatively in literatures. 

To design the electric lens devices of graphene, it is necessary to analytically study the 
focusing phenomenon of electron flow emitted from a point-like source in a graphene sheet, 
present work will fulfill the task. We will provide the analytical expressions for 2D current 
density distribution in graphene PNJ and NPNJ. Moveover, we theoretically calculated the 
electric current collecting efficiency of a drain electrode in the focus for symmetric PNJ and 
NPNJ, which play the key role in the functions of the devices. 

We start with the electronic band structure of mono-layered graphene. In a graphene, 
carbon atoms are arranged in a honeycomb lattice, each atom is connected with three nearest 
neighboring atoms by covalent bonds. The Pz electron of carbon atoms can form delocalized 
7r(7r*) band, which is valence (conduction) band of graphene respectively. The two bands 
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touch each other at six corners in the hexagonal Brillouin zone. For low-energy excitation 
near these corners, dispersion relation of electron is linear as Ec(v){k) = ^hvp \k\, which is 
closely similar to the energy spectrum of 2D massless Dirac fermionsii"— . The group velocity 
of electron is Vc(v) = dE {k^d{hk^ = ±vpk/k (" +" for conduction band electron, and " — " 
for valence band electron). The negative refractory index phenomenon is a consequence of 
different signs in the expressions of Vc{v) on two sides of the PNJ and the conservation of 
momentum of wave in the direction perpendicular to the direction that the potential varies. 

Ignoring the inter-valley scattering and the degree of freedom for spin, the low energy 
electron in graphene can be described by a two-component spinor— . The effective Hamil- 
tonian is H = hvpa ■ V, with vp ~ c/300 = 10^ cm • s^^, a = (cr^, cXy), and ay are Pauli 
matrices. We consider a circular source electrode with the center being located at (—a, 0). 
The radius of electrode should be much smaller than the size of whole graphene sheet, but 
still satisfying kcR = 2ttR/Xf ^ 1, with Xp being the wavelength of electron wave near 
the source electrode. The wave function of zero energy electrons near the electrode can be 
explicitly described by a two-component spinor, 

I zifj^^ (A;,ro) e*'^^^'^) / V 7rA;,ro \e^H^,y) i 
which is one of a set of cylindrical function solutions of Dirac equation in left N region, 

{-iaA - ioydy - kc) tpi (x, y) = 0, (2) 

where, A is a constant to be determined by certain boundary condition, and Hn\z) is 
the Hankel function of the first kind. The iaiicj) = y/{x + 1), tq = \J {x + of' + y"^, (see 
Fig. [T]). Here we have used the lowest order asymptotic expansion of Hn\z), which implies 
far field condition kcVo 3> 1™ has been used^. It is easy to verify that only the current 
density of this type of wave function is isotropic near the edge of the source electrode, jin = 
evp'ipiaipi = lev p A? I [j^kcK) (cos0, sin0), while all other cylindrical function solutions give 
3in ~ (cos 720, sinn0), n > 2. Given the input current Jq, the continuity of electric current 
at the boundary of the source electrode implies Jq = 27ri?|jj„|, which determines constant 
A = ^k,hl{%evp). 

The penetration problem of an incident plane wave of electron at interfaces of PNJ and 
NPNJ in graphene has been already studied in literature^!^. To use these fruits in present 
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case of cylindrical incident wave, the key point is to find the suitable expansion of cylindrical 
wave in a set of plane wave basis, i.e., describe the cylindrical wave as the superposition of a 
series of plane waves. Then we will use the results for refractory (for PNJ) and transmission 
(for NPNJ) of each plane wave components of an incident wave, then add them together to 
obtain the whole refractive (transmitted) electron wave. 

The expansion of an incident cylindrical electron wave as given in Eq. can be rewritten 



asil. 



^l(ri) = -/ gifccacose, \ik4x COS 9.+y sine.) ^Q^^ ^3^ 



In general, the contour of the integral is shown in Fig. [2]^a). Here, k = {kc cos 6^, kcSinOc) is 
the wave vector of incident plane electron wave component with incident angle 6c- 

At first, we study the refractory of PNJ, which is shown in the sketch of Fig. [H We know 
a plane wave with the angle of incidence 6c has a form^^ 

/ _ ^ I "'^ 1 „ikc{x cos dc+ysindc) 

y plane— in \ .„ \ ^ i 

which refiects and refracts at the interface of PNJ, with the refiected wave and refractive 
wave being 



and 



„/, „A _ { ^ \ -iky (x cos e^+y sin e^) 

Y plane— ret r\-'^ 1 y ) /t: \ .„ \ ^ i 

V2 \ e'^'^ 



respectively. Here, three wave functions ippiane-in{x,y), 'ippiane-refi{x,y) and ippiane-refr{x,y) 
are the plane wave solutions of Dirac equations 

{-i(^xdx - icTydy + kc) tpix, y) = 0, x < 

(4) 

{-iaA - icTydy - k^) ip{x, y) = 0, x > 0. 

The factor e*'" in refiected wave function ippiane-refi{x,y) is related to Berry phased. 
By the continuity of wave function in N region and P region, 

"^plane—in 

(0, y) + rippiane-refr{0, y) = tlpplane-refl{0, y) , 
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the coefficient of refraction can be obtained easily as t = 2 cos 9c/ (e*^" +e~*^'=), and the 
coefficient of reflection is therefore r = 1 — t. 

The refractive wave function of cylindrical wave is therefore, 
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^~ik^ COS e,X-iK My ^Q^^ ^ ^ Q_ 



Momentum conservation in y direction requires kcSmOc = —kySinOy, which implies the 
relation between the angle of refraction 6^ and the angle of incidence 6^, which is similar to 
the Snell's law in optical refractory^, 



sin 9c 
sin 9„ 



n, 



(6) 



For simplification, we study the refractory for symmetric PNJ at ffist, then study the 
asymmetric PNJ with n ^ 1. For convenience, we shift the integral contour as shown in 
Fig. |2](b) for a symmetric case, i.e., = kc, n = —1, and 9^ = —9c, the refractive wave at 
P region has a simple form. 



refr 



A 



1 / e" 

d9c —p^ cos 9c 

. V2 



A 



-i{kcri--K/i) 



^ikcri cos{dc — (n—a)) 



COS a, a E [— |, |] 



(7) 



(8) 



'* nkcri 



COS a, a e [f , 



Here, tana = y/{x — a), and ri = \/ {x — a)^ + y"^. We have used steepest descent method 
with the inclusion of only one leading term^^, which implies the approximate refractive wave 
function in Eq. (|8]) is valid for k^ri ^ 1. We will explain later that present approximation 
is reasonable. The current density j = evFiplcfr^i^refr is therefore 



ji(ri,a) (k^R) 



cos^ a (cos a, sina) , a G [— 7r/2, 7r/2), 



(9) 



— cos^ a (cos a, sin a) , a G [7r/2, 3n/2). 

The intensity of current density around but not too close to the focus (a, 0) has an expression 
ji oc cos^ «/ti = (x — a)^ / [(x — a)^ + y^]^^^, as shown in Fig. [3t^a), which is in agreement 
with the corresponding result in Ref. jsj. 



If a detecting circular electrode with radius R is placed at (a, 0), the maximal electric 
current can be collected by it can be obtained as, 

|.37r/2 ^ J 

Ii = R j{R,a) ■ (cos a, sin a) = — . (10) 

A/2 4 

Thus only 1/4 of the total current from source electrode can be collected by the detecting 
drain electrode. A device of graphene electronic lens should be sophisticatedly designed for 
high efficiency. 

In general, an asymmetric PNJ leads to a refractory index n ^ —1, the refractive wave 
function does not have a simple form. The integral in Eq. ([5]) can still be calculated by 
steepest descent method. The results for n = —1.2 are shown in Fig. ISTb). As discussed in 

nn 

Refs. |3IJ18I|. electron flow transporting across an asymmetric PNJ forms caustics, with the 
cusp located at {\n\a, 0). Our results shown in Fig. 3(b) are in good agreement with those 
inRef. [b|. 

In the original prediction of optical Veselago lens, the electromagnetic wave from a point- 
like source focuses on a point after transmitting across negative refractory material, which 
is in close analogy with the focusing of electron flow at the second N region in the graphene 
NPNJ, as shown in Fig. HI 

Using similar method for refractive wave of graphene PNJ, we flrstly flnd the each plane 
electron wave component's of the transmission wave, then "add" all the components together 
(if fact, by an integral). The plane wave's transmission across the P region between two N 
regions can be solved by standard method for electron wave transmitting across a square 
energy barrier— 

We separate the 2D space in graphene NPNJ into three regions x<0, 0<x<d and 
X > d, labeled by I, II, and III, respectively. The Dirac equations for three regions of NPNJ 
are, 

[—icTxdx — icTydy + kc) V^(x, y) = 0, X < or X > d ^ ^ 

[-icTxdx - ic^ydy - K) ij{x, y) = 0, < d. 
The plane wave function at each region is the solution of Eq. (fTTj) . which can be written 
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as 



1 1 


( 1 
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a 














1 gi^c 



^ifec cos dcX+ikc sin 0c2/ 



+ 



^—ikv cos 9vX—ikv sind^y _|_ 



V2 i -e-^o. 
b ( 1 



-ifec cos (^ca:+ifcc sin 9cj/ 



, X < 0, 



^ife^ cos^^a;— jfcc sinS^y 



^ifcc cos 0ca;+ifec sin $c2/ x "> d 



, < X < d, 



The coefficients of transmission and reflection, as well as coefficients a and b, can be 
obtained by the continuity of wave functions at two interfaces x — and x — d, 

V'i(0, y) = -011(0, y), ijjuid, y) = ?/)■ 

For symmetric NPNJ, n — —1, 9y — —9c, and A;^, = kc, the coefficient of transmission 

cos^ 9r 



t 



cos 9re~ 



g2i feed cos 9c sin^ ^ 

Thus, the transmission wave at the second N region can be written as, 

1 



r-n/2—ioo 

'4^transmissioni,-^ y y^ — / COS 9c G 

^ J-TT/2+ioo \ e" 



^kcvM^Q^^ (12) 



which can be rewritten as. 



transmission 



{x,y) 



j\ r2~^°° g2*yc _|_ g-2«e;c _|_ 2 / 1 



oo 



Here, r2 = \J[^ ~ (2c? — a)]^ + y^, and tan^ = — {2d — a)]. Using the steepest descent 
method, we have. 



'0traiismission(2^, ?/) ~ ^ 



■KkcTi 



Trkcr2 



COS^ l3 



COS^/3 



3j(fccr2-7r/4) 



/5 e [-1, f] 



(14) 



il3 



TT StT 
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Here, ai = 2d — a. At the second N region (Region III), the electric current density around 
the point (ai, 0) is 

J(r2, /3) = ± \""''^ (cos/3, sin/3), fc.ra » 1, , (15) 

Here, "+" for /3 G [-7r/2, 7r/2], and "-" for /3 G (7r/2, 37r/2], similar to those used in Eq. 
Thus around the focus but not too close to it, the intensity of current density is 

32 = \j{x,y)\ 



2- [ix-a,r + ff'' 

as shown in Fig. |5l The maximal current can be collected by a drain electrode at focus 
(ai, 0) is /q /Jjg''^ cos"^ /3/(27r) d/3 = 3/16 Jq, which is nearly one-fourth of the total current 
radiated from source electrode. 

In above analytical expressions, we keep only one leading term in the results obtained 
by the deepest descent method, which implies that we consider the wave at the positions 
ri ^ Xp (NPJ) and r2 ^ Xf (NPNJ). We know the low-energy effective Hamiltonian 
of graphene is valid for energy less than 1 eV^^, thus is at most several nanometers. 
In experiments of graphene-based devices, the source/ drain (detecting) electrode's size is 
roughly tens of nanometers, thus this far field condition k^ri {kcr2) ^ 1 is quite reasonable. 
In fact, due to the short wavelength of electron wave, high resolution for one wavelength or 
sub- wavelength is less meaningful than that for photons^. 

In summary, we have theoretically reproduced the electric lens phenomena in graphene 
PNJ and NPNJ, and have obtained the analytically expressions for the current density 
distribution in PNJ and NPNJ. The key idea is to expand the incident cylindrical wave in a 
series of plane wave basis, then find the refractive wave or transmitted wave for each incident 
plane wave component, adding them together, get the total refractive wave (in PNJ) and 
transmitted wave (in NPNJ), as shown in Eqs. fl9|)f lT5|) . Fig. [3l and Fig. [51 The analytical 
results are in good according with available numerical results in Ref.-. We firstly obtained 
the maximal possible current which can be collected by drain electrode in symmetric PNJ 
and NPNJ, which are 1/4 and 3/16 respectively. These data are important in designing 
devices, such as logical gates and interconnects based on graphene Veselago lens. 

A recent progress on the fiat-lens focusing of electrons on the surface of topological 
insulator suggested a high efficient electric lens effect based on topological insulator—. In this 
electric lens, only conduction band electrons are used, which avoids high interface resistance 
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in graphene PNJ. We are now working on the possibility of applying our theory into this 
interesting issue. 
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FIG. 1: The schematic diagram for electric lens in symmetric graphene PNJ. The electric current 
source is located at {—a, 0). The electron flow radiated from current source focuses at P region. 
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FIG. 2: (a) The general integral contour of Eq. [5] in complex plane of complex variable 9c— and 
(b) The special integral contour for symmetric PNJ with refractory index n = — 1. 
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FIG. 3: Intensity distribution of current density in P region of a graphene PNJ {x > 0), x and y 
coordinates are in unit of a. (a) n = —1, focus is located at (a,0), (b) n = —1.2, refracted wave 
forms caustic, with the cusp located at (1.2a, 0). 
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FIG. 4: Schematic diagram for symmetric electric lens in graphene NPNJ. 



12 




-U.D 



-1.0 



-0.5 



0.0 



0.5 



1.0 



1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



FIG. 5: Intensity distribution of current density in second N region (x > d) of a graphene 
The X and y coordinates are measured from the focus located at {2d — a, 0), in unit of a. 
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